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Abstract 

The filling of tanks with cryogens In the low- 
gravity environment of space poses many technical 
challenges. Chief among these Is the Inability to 
vent only vapor from the tank as the filling pro- 
ceeds. As a potential solution to this problem, 
the NASA Lewis Research Center (NASA Lewis) Is 
researching a technique known as No-Vent Fill. 

This technology potentially has broad application. 
The focus of this paper Is the fueling of space 
based Orbital Transfer Vehicles. 

This paper describes the fundamental thermo- 
dynamics of the No-Vent Fill process to develop an 
analytical model of No-Vent Fill. The model is 
then used to conduct a parametric investigation of 
the key parameters; initial tank wall temperature, 
liquid-vapor Interface heat transfer rate, liquid 
Inflow rate and inflowing liquid temperatures. 
Liquid inflowing temperature and the liquid-vapor 
Interface heat transfer rate seem to be the most 
significant since they Influence the entire fill 
process. The initial tank wall temperature must be 
sufficiently low to prevent a rapid pressure rise 
during the Initial liquid flashing stage, but then 
becomes less significant. 

Nomenclature 

A area 

Cv specific heat at constant volume 

h enthalpy 

h convective heat transfer coefficient 

k thermal Conductivity 

L characteristic length 

M mass 

m mass flow rate 

q heat flux 

Q heat Input 

T temperature 

t time 

U total Internal energy 

u Internal energy per unit mass 

V volume 

N work 

N rate of work 
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a 

thermal diffusivity 

P 

density 

Subscrl 

pts: 

cond 

condensate 

fgas 

flashing gas 

f 1 1 q 

flashing liquid 

gas 

ullage gas 

In 

inlet 

inf 

Interface 

ig 

between liquid and gas 

1 1 q 

bulk liquid 

sat 

saturation 

sgas 

saturated gas 

s 1 1 q 

saturated liquid 

tank 

tank 

wal 1 

tank wal 1 


Introduction 

Many fluid management Issues face the designer of 
systems for filling tanks with cryogens on-orbit. 
The most Important of these are the large rate of 
generation of vapor from the residual energy stored 
In cryogenic tank walls, the uncertainty of liquid 
and vapor distributions in a tank In low gravity, 
and the need to keep maximum tank pressure low to 
reduce tank mass. During a normal gravity fill a 
top vent is kept open to vent the vapor generated 
during the fill process thereby maintaining a low 
tank pressure. If the same approach is used in low 
gravity the ullage gas may never vent. Instead of 
venting vapor, large amounts of liquid may be 
dumped overboard. If liquid Is vented from one 
side of a nonpropul si ve vent and vapor from the 
other side, the spacecraft may tumble out of con- 
trol. To prevent the loss of liquid, the space- 
craft may be placed in an artificial gravity field 
by continuous thruster firing to position the 
ullage at a vent opening, but this may require 
dedicated thrusters and additional propellent. 

The problems of low-g fluid management and the 
thruster technique for liquid positioning have 
been described previously by Lacovlc. 1 For the 
fill operation, maintaining control between the 
tanker spacecraft and the spacecraft being refueled 
while thrusting may be difficult due to potentially 
large shifts In the center of gravity. In some 
Instances, such as a depot based at a space sta- 
tion, thrusting may be Impractical due to large 
system size. The No-Vent Fill technique, first 
analyzed in detail by Merino, Is a methodology 
whereby these problems are eliminated. 2 



The No-Vent Fill methodology is as follows: 

The tank wall is prechilled to a temperature suffi- 
cient to remove most of the thermal energy from it. 
The tank is evacuated by venting to space, and liq- 
uid is then injected through spray nozzles and/or 
mixing jets to promote liquid-vapor heat exchange. 
The tank can be filled without venting as long as 
the prechill Is sufficient and the liquid-vapor 
heat exchange keeps the liquid and vapor close to 
thermodynamic equilibrium. 

In the main body of this paper, a model for 
the initial flashing stage (when tank pressure Is 
less than the liquid saturation pressure) is formu- 
lated. A model for the later stages, which pre- 
serves the temperature difference between the bulk 
liquid and the vapor, is also formulated. These 
two models are combined to provide a transient 
model of the No-Vent Fill process. A computer pro- 
gram was developed to solve the model equations, 
and then used to conduct a parametric study of the 
No-Vent Fill process. Parameters Investigated 
include: initial tank wall temperature, interface 

heat transfer, liquid Inflow rate, and inflowing 
liquid temperature. 

No-Vent Fill Theory 


General 


Following the tank chilldown and evacuation, 
liquid inflow is started. Liquid is Injected 
through spray nozzles and/or mixing jets to promote 
liquid-vapor heat exchange. The initial Inflow 
partially flashes while entering, and vaporizes on 
striking the tank wall. This removes the remaining 
residual energy from the tank walls and raises the 
tank pressure to near saturation conditions. The 
continuing Inflow of liquid either condenses or 
compresses the vapor allowing more liquid to enter. 
For a period of time condensation dominates and 
tank pressure remains fairly constant. Eventually 
vapor compression takes over. The tank pressure 
rises rapidly to the limit pressure, and the fill 
process must be stopped. It is the objective of 
No-Vent Fill technique to postpone the compression 
phase as long as possible. 


Liquid Flashing Stage 

The initial flashing stage is modeled as fol- 
lows. The high heat transfer rates found in boil- 
ing make it unlikely that much liquid will remain 
in the tank until the wall is chilled to near the 
incoming liquid temperature. When the tank pres- 
sure is less than the saturation pressure of the 
incoming liquid, the liquid Inflow will flash 
resulting in a gaseous inflow which enters the 
ullage vapor directly and a liquid stream which Is 
cooled to the saturated state corresponding to the 
current tank pressure. This liquid stream then 
strikes the tank wall and is assumed to be quickly 
vaporized to a saturated vapor. A comparison 
between the heat transfer rates of film boiling of 
the liquid stream and gaseous conduction Into the 
ullage gas (the comparison is included in an appen- 
dix to this report) indicates that the effects of 
wall-gas heat transfer are secondary as long as 
conduction is the predominate heat transfer mechan- 
ism. For simplicity the effects of wall-gas heat 
transfer will be assumed negligible. The formula- 
tion presented here only accounts for the tank wall 
energy removed by the vaporizing liquid stream. 

The equations for this stage are: 

Inflow energy balance at the liquid inlet 

m.h, — mph + mp»*h<i (1) 

in in fgas sgas f 1 1 q sliq 

Initial flashing mass balance 

"in ' "fgas* "fllq (2> 

Wall energy balance 

d(C v T ) . 

M wal 1 dt = "fl1q <h sgas " h sl1q ) (3> 
Gas mass balance 


dM 

gas ■ • 

dt ' "fgas + "fllq 


Gas energy balance 


(4) 


To understand the No-Vent Fill concept it is 
helpful to examine the thermodynamics of the Fill 
process. It can be shown by an equilibrium analy- 
sis of the thermodynamics that, with a sufficiently 
low temperature of the entering liquid, the end 
state of a full tank should be at an acceptably low 
pressure. Unfortunately like all real processes, 
the fill process does not proceed at thermodynamic 
equilibrium. Major sources of nonequilibrium 
behavior are as follows: (1) energy transfer from 

the initially warm tank to the entering fluid, (2) 
heat flux Into the tank during the fill process, 

(3) flashing of the liquid during the Initial low 
pressure stages of the fill, and (4) temperature 
gradients within and between the tank, the vapor, 
and the liquid. 

A model which separates the vapor and liquid 
phases is needed to determine the effect of temper- 
ature differences and the rate of heat transfer 
between the two phases. Since the initial flash- 
ing stage of the flow is controlled by different 
phenomena than the condensation and compression 
stages, the model will be divided into two parts. 
For simplicity, heat flux into the tank during the 
fill process will be assumed to be negligible. 


d(M na U a ) 

3 dF 9 = ( "fgas * m fl1q )h sgas <5) 

Rearranging Eq. (1) and substituting Eq. (2) into 
Eq. (1) 

"ln h 1n = ( "in " "fl1q )h sgas = "s11q h fliq (6) 
or 

m* (h. — h ) = nip , , ( h * , — h ) (7) 

in in gas fllq sliq sgas 

Substituting Eq. (7) Into (3) 

d(C v T> . 

" M wall di = "ln <h sgas " h 1n> <8) 


Substituting Eq. (2) into (4) 
dM 

911 _ « 


dt 


( 9 ) 


and into Eq. (5) 
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d(M 


qas u gas ) 

dt 


" ™in h sgas 


( 10 ) 


Expanding Eq. (5) and using the results of (9) and 
( 10 ): 


du 

A = m. Oi - u rt . 

gas dt In sgas gas) 


( 11 ) 


For the times when the tank pressure is above 
the inflowing liquid saturation pressure, Eqs. (8), 
(9), and (11) can be written directly as the mass 
and energy equations for the Fill process. Once 
the liquid and wall temperatures are equal the 
residual wall energy is negligible and the model 
can be switched to the condensation and compression 
stage equations. 


this document. The main effects of the gas Is the 
mass transport of energy already accounted for in 
Eq. (17). qj n f is then a function of liquid side 
heat transfer. Writing this as a convective 
equation: 


qinf - hA lnf< T 1nf - T 1 1 
It can be shown by kinetic theory that: 

Tlnf ■ T sat • 

Combining Eqs. (17) to (19). 

^ A 1nf (T sat " T 1iq ) 
cond = (h ga$ - h liq ) 


(18) 

(19) 


(20) 


Vapor Condensation and Compression Stage 

To study the relation between gas and liquid 
the fluid is broken into three nodes: gas, liquid 

and interface. A mass and energy balance can be 
written for each node. 


The work on the gas can be calculated from basic 
thermodynamics for a control volume: 

W = PdV (21) 

If the pressure rise is not too great: 


Gas mass balance: 


dM nac 

— - = _m 

dt cond 


( 12 ) 


Gas energy balance: 

d(M qas u qas) H du qas 
dt s M gas dt" 


+ m cond u gas 
^cond^gas + ^lg 


Liquid mass balance: 
dM 


115 


dt 


“in 


+ m 


cond 


(13) 


(14) 


Liquid energy balance: 


du 


M 


115 


dM 


liq dt 


+ u. 


115 


dt 


+ + ril 1n h 1n + m 


w 


ig 


cond h iq 
(15) 


The Interface node has some unique features. 

It Is assumed to be a infinitesimally thin layer 
between the liquid and gas. Since It Is Infinites- 
imally thin, it can store neither mass nor energy. 
This makes the equations for the Interface as 
follows: 


Interface mass balance: 
dM 


inf 


dt 


= 0 


Interface energy balance: 


q 1 nf + m cond h 11q 


"'cond h gas 


Some additional equations are needed to solve 
this set of equations. qj n f is determined by a 
convective heat transfer equation. Heat transfer 
from the gas side can be neglected because gas 
thermal conductivity Is much smaller than liquid 
thermal conductivity for the fluids of Interest In 


dV 

W - P -9 as - 
w lg gas dt 

In a tank of constant volume: 


V r = V,. = constant 

gas liq 


From Eq. (23): 


dV na( . dV,. 
gas liq 

dt = dt 


From Eq. (24) 
dV 


Ha 


dt 


p llq 


(m 


1 n 


^cond^ 


So: 


W, = (lii, + m 
lg in cond 


( 22 ) 


(23) 


(24) 


(25) 


(26) 


(16) 


(17) 


Equations (12) to (17), (20), and (26) along 
with the thermodynamic equations of state form a 
fairly complete model of the No-Vent Fill process. 
The chief source of thermodynamic inefficiency Is 
the work term W] g . The maj<?r independent vari- 
ables in this analysis are m^ n , h\ n , T wa n , and 

Vtank- 

Anal ys 1 s 

The equations for the two stage model are suf- 
ficiently complex that a closed form solution Is 
not possible. A computer algorithm called NVFILL 
has been written to solve the equations using a 
finite difference approximation. The variables 
min, hi n , h, and Vtank are inputs to the code 
and held constant throughout a run. T wa n Is set 
at the start. A^ n f Is defined as the surface 


area of a sphere whose volume Is that of 
the gas volume which changes with time. 


‘'gas* 


During the initial flashing stage, the compu- 
tational algorithm adds the mass added to the 
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ullage during the timestep, calculates a new 
ullage density Internal energy and enthalpy, and 
then uses the density and enthalpy to determine a 
new ullage pressure. In the condensation- 
compression stage the problem is solved by a two- 
step calculation procedure. For the first step, 
the ullage is held at constant pressure while a 
interfacial mass transfer rate is calculated. For 
the second step, a new liquid volume is calculated 
from the liquid inflow, interfacial mass transfer 
and current bulk liquid conditions. The ullage 
vapor is then compressed ad i aba t 1 ca 11 y to fill the 
remaining tank volume. The time step for both 
these processes combined is set to 10 -3 hr which 
is sufficiently small compared to the process rates 
for most cases of interest to insure a good approx- 
imation of the continuous liquid inflow, interfa- 
cial mass transfer and vapor compression processes. 
The computational algorithm terminates when the 
tank is 95 percent full with liquid or the tank 
pressure exceeds 60 psia. 

Results and Discussion 

Many differing tank sizes and geometries exist 
in the literature for Orbit Transfer Vehicle (OTV) 
designs. To provide a comparison with previous 
work, the OTV tank parameters employed by DeFelice 
and Aydelott are used. 3 These OTV parameters are 
summarized in Table 1. The baseline OTV parameters 
in Table 1 are representative of an unmanned space 
based OTV and are expected to be the sizes used in 
an initial operational space based OTV. Only the 
hydrogen tankage is examined since previous work 
indicates that the fill transients in the hydrogen 
tankage are the most severe. 2 

The effect of initial tank temperature on the 
tank pressure, as filling proceeds, Is shown in 
Fig. 1. The 218.3 K <393 °R> temperature was the 
suggested Initial tank temperature from Ref. 3. 

It produces a large initial pressure transient and 
causes the tank to fill at a much higher pressure 
than the other initial tank temperatures investi- 
gated. The 152.8 K (275 °R) initial condition was 
selected as a wall temperature which will raise the 
pressure in the tank just above the stagnation 
pressure of the incoming liquid at the end of the 
vapor flashing stage. A temperature of 55.6 K 
(100 °R) was selected as a likely minimum tempera- 
ture, below which little effect on tank pressure 
results because all the tank energy has been essen- 
tially removed. It can be seen from Fig. 1 that 
even though the 55.6 K (100 °R) initial tank tem- 
perature transitions from the wall chilldown to the 
compression/condensation stage at about 130 kPa 
(15 psia) bulk boiling in the liquid raises the 
tank fill pressure characteristics to near the 
152.8 K (275 °R) line for the majority of the fill 
process . 

Reference 3 selected the Initial tank tempera- 
ture based on a criteria which yielded the desired 
fill level at the end of the fill using a thermo- 
dynamic equilibrium analysis to predict the end 
state. The initial pressure transient It produces 
was not accounted for in the Ref. 3 analysis and 
may be undesirable because of the control problems 
it poses for tanks with strict maximum pressure 
requirements. The stagnation pressure requirement 
used to select the intermediate Initial wall tem- 
perature is probably a better criteria for initial 
wall temperature since it minimizes the initial 
pressure transient. Consequently, the 152.8 K 


(275 °R) wall temperature was used for the rest of 
the No-Vent Fill studies. 

Previous work on fluid mixing give wide ranges 
of surface heat transfer coefficients. 4-6 To cover 
the range of expected values, heat transfer coeffi- 
cients spanning the range of two decades from 5.68 
to 568 W/m 2 K (1 to 100 Btu/ft 2 hr °R) were para- 
metrically evaluated. As can be seen in Fig. 2 
(data from Fig. 1 assumed 56.8 W/m 2 K 
[10 Btu/ft 2 hr °R] ) , heat transfer coefficient has 
a significant effect on the tank pressure response 
during the No-Vent Fill process. For the smallest 
heat transfer coefficient, tank pressure exceeds 
414 kPa (60 psia) prior to completion of the fill 
process. This is the only analysis run which did 
not terminate at 95 percent full. The highest heat 
transfer coefficient produces hardly any pressure 
rise beyond the initial flashing stage. 

The 56.8 W/m 2 K (10 Btu/ft 2 hr °R) heat trans- 
fer coefficient was selected as a estimate for No- 
Vent Fill calculations for the rest of the study. 
Heat Transfer rates above this are likely to be 
achieved only by using forced convection to promote 
liquid-vapor heat transfer. 

Figure 3 shows the tank pressure response for 
No-Vent Fills at liquid inflow rates at 453.5 kg/hr 
(1000 Ibm/hr), 907 kg/hr (2000 lbm/hr, Table 1 
rate) and 1814 kg/hr (4000 lbm/hr). Although the 
curves are fairly close for most of the fill proc- 
ess there Is a significant difference in the final 
pressure with the highest pressure corresponding 
to the highest flow rate. 

Figure 4 shows fills at various inflowing liq- 
uid temperatures corresponding to saturation pres- 
sures of 34.5 kPa (5 psia), 103 kPa (15 psia) and 
172 kPa (25 psia). Inflowing liquid temperature 
shows a very strong effect on the No-Vent Fill 
process. This is to be expected given the depend- 
ence of the condensation rate on temperature dif- 
ference. For real spacecraft, subatmospheri c 
subcooling (lowering the liquid temperature below 
the saturation temperature at atmospheric pressure) 
may be difficult to achieve so 103 kPa (15 psia) 
was used as the baseline for the other parametric 
runs. The pronounced effect of subcooling indi- 
cates that further investigation of heat exchangers 
to cool the liquid Inflow is probably worthwhile. 

Concluding Remarks 

Several areas for model Improvement can be 
suggested as a result of this study. Most evident 
Is the need for a more accurate correlation relat- 
ing the Interfacial heat transfer rate to the Inde- 
pendent parameters of No-Vent Fill. Two different 
correlations are under investigation currently but 
are dependent on liquid injection technique; one 
for droplet sprays through the ullage, and a sec- 
ond for turbulent jet mixing of the bulk liquid. 

In most real transfer systems the Inflow rate will 
not be constant but will be controlled by the pres- 
sure difference between the liquid supply and the 
tank being filled. A model which links the flow 
rate to the tank pressure variation is being deve- 
loped. The effects of wall-vapor and vapor-liquid 
heat transfer on the initial vapor flashing stage 
which were neglected in the analysis of this paper 
are under continued investigation. Several pro- 
posed injection techniques for No-Vent Fill show 
the potential for promoting forced convection in 
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the vapor thus raising the wall vapor heat transfer 
rate much closer to the liquid boiling heat trans- 
fer rate. The exact heat transfer correlations for 
these conditions are still being determined. 

The results of these parametric investigations 
help define areas of research where the No-Vent 
Fill process should be studied. Several ground 
test experiments are being prepared at NASA Lewis 
to study the controlling phenomena. One experimen- 
tal apparatus will employ liquid Nitrogen as the 
test fluid and provide controlled flows of liquid 
and gas which can be used to study heat transfer 
and condensation on the liquid interface. A small 
portable experimental test rig capable of being 
operated in a remote area with LH 2 is also being 
fabricated. This test rig will be used to study 
the effect of thermal subcooling and injection 
technique on the No-Vent Fill process. In addi- 
tion, NASA Lewis is reactivating a large cryogenic 
vacuum chamber to allow No-Vent Fill testing with 
LH 2 on large scale tankage. Initial testing for 
this facility will investigate the effects of ini- 
tial tank temperatures and in flow rate, as well as 
provide proof of concept testing for the No-Vent 
Fill technique. It is hoped that the combined 
results from the portable and large scale tests 
will provide the data necessary to determine inter- 
facial heat and mass transfer rates and verify the 
No-Vent Fill analytical models. 

Although ground testing can be used to inves- 
tigate key phenomena, it is only in the low- 
gravity environment where the fluid dynamics and 
heat transfer mechanisms combine to give a true 
simulation of the Low-Gravity No-Vent Fill process. 
To investigate the effect of a Low-Gravity environ- 
ments, several approaches are either under develop- 
ment or being planned. NASA Lewis has had good 
success with Computational Fluid Dynamic codes to 
model low-gravity flows. An effort is underway to 
adapt these codes to the No-Vent Fill problems. 

NASA Lewis is undertaking studies for the design 
of a free flying spacecraft called the Cryogenic 
Orbiting Liquid Depot Satellite (COLD-SAT) to 
study fluid management issues with liquid hydrogen 
in space. The study of No-Vent Fill is one of the 
primary experiments to be conducted on the 
COLD-SAT. Although a LH 2 experiment such as 
COLD-SAT is required to understand the technology 
completely, NASA Lewis is also exploring small 
scale low-gravity experiments which can provide 
precursory data for COLD-SAT and allow for a more 
focused COLD-SAT design. 

Appendix 

To estimate the heat transferred from the tank 
wall to the ullage vapor by means of solid conduc- 
tion, the vapor is approximated as a semi-infinite 
slab of gaseous hydrogen. Schneider gives as the 
solution for the total heat input for a semi- 
infinite slab. 7 



Gaseous Hydrogen at 101.3 kPa (15 psla) has the 
following property values: 8 

T gas = 20.268 K (36.6 °R> (A2) 

k = 0.01694 W/K m (0.00979 Btu/ft 2 hr °R> (A3) 


a = 0.00375 m 2 /hr (0.04037 ft 2 /hr) (A4) 
For a typical fi 1 1 

T wa n = 152.8 K (275 °R) (A5) 

t = 0.055 hr (A6) 

(for the duration of the flashing stage.) 

Substituting these values into Eq. (Al) 

2 = 9.69 N hr/m 2 (3.07 Btu/ft 2 ) (A7) 

Dividing by the time of the flashing stage to 
obtain an average heat flux 

3 = 176 W/m 2 (55.8 Btu/ft 2 hr) (A8) 


Brentari and Smith give as the minimum film boil- 
ing heat flux for liquid hydrogen a value of 
approximately: 9 

3 = 8.5 x 10 3 W/m 2 (2.69xl0 3 Btu/ft 2 hr) (A9) 

The gaseous conduction heat flux is therefore only 
2.1 percent of the minimum heat flux to be found 
in the film boiling process. 
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Table 1. - Prototype OTV LH 2 Tankage Characteristics 


Mass, kg (lbm) 256 (565) 

Volume, m 3 (ft 3 ) 43 (1500) 

Mass to volume ratio, kg/m 3 (lbm/ft 3 ) .... 6 (0.37) 

Inflow rate, kg/hr (lbm/hr) 970 (2000) 


INITIAL 

WALL 



TIME FROM START OF FILL, HR 

FIGURE 1. - EFFECT OF INITIAL WALL TEMPERATURE ON THE NO- 
VENT FILL PROCESS. PROTOTYPE OTV LH2 TANKAGE; TANK VOL- 
UME 1500 cu FT; MASS TO VOLUME, 0.3 LBM/cu ft; INFLOW, 
2000 lbm/hr; INFLOW TEMPERATURE, 36. G R; HEAT TRANSFER 
COEFFICIENT, 10 Bty/sq ft-hr-R. 
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HEAT 

TRANSFER 

COEFFICIENT, 



FIGURE 2. - EFFECT OF INTERFACIAL HEAT TRANSFER CO- 
EFFICIENT ON THE NO-VENT FILL PROCESS. PROTOTYPE 
OTV LH2 TANKAGE; TANK VOLUME 1500 CU FT; MASS TO 
VOLUME 0.3 LBM/CU FT; INFLOW, 2000 LBM/HR; INFLOW 
TEMPERATURE, 36.6 R; INITIAL TANK WALL TEMPERATURE. 
275 R. 


LIQUID INFLOW 



TANK LIQUID FILL, PERCENT 


FIGURE 3. - EFFECT OF LIQUID INFLOW RATE ON THE NO-VENT 
FILL PROCESS. PROTOTYPE OTV LH2 TANKAGE; TANK VOLUME 
1500 CU FT; MASS TO VOLUME, 0.3 LBM/CU FT; INITIAL WALL 
TEMPERATURE, 275 R; INFLOW TEMPERATURE, 36.6 R; HEAT 
TRANSFER COEFFICIENT, 10 Btu/sq ft-hr-R. 
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INFLOW 

TEMPERATURE, 



FIGURE H. - EFFECT OF INFLOWING LIQUID TEMPERATURE ON THE 
NO- VENT FILL PROCESS. PROTOTYPE OTV TANKAGE; TANK 
VOLUME 1500 cu FT; MASS TO VOLUME, 0.3 LBM/cu FT; 
INFLOW, 2000 LBM/HR; INITIAL WALL TEMPERATURE, 275 R; 
HEAT TRANSFER COEFFICIENT, 10 Btu/sq FT-HR-R. 
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